Heavy Quark Energy Loss in Nuclear Medium 
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Multiple scattering, modified fragmentation functions and radiative energy loss of a heavy quark 
propagating in a nuclear medium are investigated in perturbative QCD. Because of the quark mass 
dependence of the gluon formation time, the medium size dependence of heavy quark energy loss 
is found to change from a linear to a quadratic form when the initial energy and momentum scale 
are increased relative to the quark mass. The radiative energy loss is also significantly suppressed 
relative to a light quark due to the suppression of collinear gluon emission by a heavy quark. 



An energetic parton propagating in a dense medium 
suffers a large amount of energy loss due to multiple 
scattering and induced gluon bremsstrahlung [1]. In a 
static medium, the total energy loss of a massless parton 
(light quark or gluon) is found to have a quadratic de- 
pendence on the medium size [2-6] due to non-Abelian 
Landau-Pomeranchuk-Migdal (LPM) interference effect. 
In an expanding medium, the total energy loss can be 
cast into a line integral weighted with local gluon density 
along the parton propagation path [7-9]. Therefore, the 
measurement of parton energy loss can be used to study 
properties of the medium similar to the technique of com- 
puted tomography. Recent experimental measurements 
[10,11] of centrality dependence of high-py hadron sup- 
pression agree very well [12] with such a parton energy 
loss mechanism. 

Because of the large mass of the heavy quark with a 
velocity v w 1 — M 2 /2E 2 , the formation time of gluon 
radiation, r/ ~ l/(u) g M 2 /2E 2 +& r /2oj g ) is reduced rel- 
ative to a light quark. One should then expect the LPM 
effect to be significantly reduced for intermediate energy 
heavy quarks. In addition, the heavy quark mass also 
suppresses gluon radiation amplitude at angles smaller 
than the ratio of the quark mass to its energy [13] rela- 
tive to the gluon radiation off a light quark. Both mass 
effects will lead to a heavy quark energy loss different 
from a light quark propagating in a dense medium. This 
might explain why one has not observed significant heavy 
quark energy loss from the PHENIX [14] measurement 
of the single electron spectrum from charm production 
in Au + Au collisions at yfs = 130 GeV. In this Letter, 
we report a study on medium induced energy loss and 
the modified fragmentation function of a heavy quark. 
In particular, we will show how the mass effects reduce 
the total energy loss and how the medium size depen- 
dence changes from a linear dependence to a quadratic 
one when the energy of the heavy quark or the momen- 
tum scale is increased. Similar results have been reported 
in Ref. [15] during the completion of this work. 

To separate the complication of heavy quark produc- 



tion and propagation, we consider a simple process of 
charm quark production via the charge-current interac- 
tion in DIS off a large nucleus. The results can be easily 
extended to heavy quark propagation in other dense me- 
dia. The differential cross section for the semi-inclusive 
process £{L{) + A(p) — > ve(L 2 ) + H(£ H ) + X can be 
expressed as 
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Here L\ and L 2 are the four momenta of the incoming 
lepton and the outgoing neutrino, £r the observed heavy 
quark meson momentum, p = [p + , m 2 N /2p + , 0j_] is the 
momentum per nucleon in the nucleus, and s = (p+Li) 2 . 
Gp is the four-fermion coupling constant and q = L 2 — 
L\ = [— Q 2 /2q~, q~, 0j_] the momentum transfer via the 
exchange of a VF-boson. The charge-current leptonic ten- 
sor is given by L« = l/2Tr(£ l7/i (l - Ts )£ 2 (l - Ts ) 7l/ ). 
We assume Q 2 <C Myy- The semi-inclusive hadronic ten- 
sor is denned as, 



Ei H |^ = \ £ (A\ J+\X, H) (A, H\J+*\A) 



x 2nd (q + p-px - £h) 
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where J2x runs over an possible final states and J+ = 
C7 M (1 — 75)50 is the hadronic charged current. Here, 
s$ = scosOc ~ dsin^c and Oc is the Cabibbo angle. To 
the leading-twist in collinear approximation, the semi- 
inclusive cross section factorizes into the product of quark 
distribution ff e {xB + Xm), the heavy quark fragmenta- 
tion function Dq^h{zh) {%h — £h/£q) an< ^ ^ ne nar d 
partonic part H$ (k,q, M) [17]. Here, xb = Q 2 /2p + q~ 
is the Bjorken variable and xm = M 2 /2p + q~. 

Similar to the case of light quark propagation in nu- 
clear medium [6] , the generalized factorization of multiple 
scattering processes [16] will be employed. We will only 
consider double parton scattering. The leading contri- 
butions are the twist-four terms that are enhanced by 
the nuclear medium in a collinear expansion, assuming a 
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small expansion parameter a s A 1 / s /Q 2 . The evaluation 
of 23 cut diagrams are similar to the case of a light quark 
[17]. The dominant contribution comes from the central 
cut diagram, giving the semi-inclusive tensor for heavy 
quark fragmentation from double quark-gluon scattering, 
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where 



T A > c (x, x Ll M 2 ) = \l ^-dy-dy-HS 
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are twist-four quark-gluon correction functions of the nu- 
cleus. Here 



Hg = Cl (z,£ 2 T ,M 2 )(l-e- lXL P 
+ c 2 (z,£ 2 T ,M 2 )[e- l7LP+y i(l 

-ix L p + (y~ - 



-ixLP + (y 



-ix L p+(y 



-Vi 



+ e 



~ Vl '(1 — e~ iXLP 2/2 



+ c 3 {z,£ T ,M 2 )e 



-ix L p + (y -y 1 ). 



-Vi )) 



(5) 



and 
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where, 2+ = Xl + (1 — z)xm/z is the additional frac- 
tional momentum of the initial quark or gluon in the 
rescattering that is required for gluon radiation, and 
xl = £^/2p + q^ z(l — z). The contribution from gluon 
fragmentation is similar to that from quark fragmenta- 
tion with z — > 1 — z. The virtual correction can be ob- 
tained via unitarity constraint. One can recover the re- 
sults for light quark rescattering [18] by setting M = in 
the above equations. Notice that we have embedded the 
phase factors from the LPM interference in the effective 
twist-four parton matrix T A ' C (x, xl,M 2 ). 



Rewriting the sum of single and double scattering 
contributions in a factorizcd form for the semi-inclusive 
hadronic tensor, one can define a modified effective frag- 
mentation function Dq^h(zh 1 M 2 ) as 
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where Dq^h( z h, M 2 ) and -D s ^ir(z#, A* 2 ) are the leading- 
twist fragmentation functions of the heavy quark and 
gluon. The modified splitting functions are given as 



A7,- M (*)= ^^f A f(x,x L ,M 2 ) 

+ 6(l-z)AT A > c (x,t 2 T ,M 2 )] 
2-nC A a s i\ 
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(l-z) 2 M 2 ]+V c /^(x)' 

AT A g c (x,£ 2 Tl M 2 ) ee f [2T A ' c (x,x L ,m 2 )\ z=1 
Jo 1 — z 

-(l + z 2 )T q A g c (x,x L ,M 2 )] , (11) 

and A-y q ^ gq (z) = Aj q ^ qg (l — z). Here we have 
suppressed other variables in A7. Given the twist- 
four quark-gluon correction functions of the nucleus, 
T A ' C (x, l\, M 2 ), one should be able to evaluate the mod- 
ified heavy-quark fragmentation function. 

As seen from the phase factors in the effective twist- 
four matrix clement Eq. (5), the gluon formation time 
for radiation from a heavy quark is 
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which is shorter than that for gluon radiation from a light 
quark. This should have significant consequences for the 
effective modified quark fragmentation function and the 
heavy quark energy loss. 

As discussed previously [19], one can assume a factor- 
ized form of the twist-four parton matrix 



Tqg C (x, xl, M 2 ) « —f q A (x){(l e-*H**) 
x A 1. 

x [ Cl (z,4,M 2 )-c 2 (z,4,M 2 )] 
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in the limit 11 C ir < ^- Here xa = I/itinRa and 
xt = (k T )/2p + q^ z is the momentum fraction associated 
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with the jnitial intrinsic transverse momentum. The co- 
efficient C = 2Cxxfg \xt) should in principle depend on 
Q 2 . With this simplified form of twist-four matrix, one 
can then calculate the heavy quark energy loss, defined 
as the fractional energy carried by the radiated gluon, 
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where xm = (1 — z)xm/z and x^ — /i 2 /2p 
xm- Note that xl/xa = L A /rf with L A 



- z) + 



the nuclear size in the chosen frame. The LPM interfer- 
ence is clearly contained in the second term of the inte- 
grand that has a suppression factor 1 — e~ XL ^ XA . The first 
term that is proportional to 03(2, M 2 ) corresponds to 
a finite contribution in the factorization limit. We have 
neglected such a term in the study of light quark propa- 
gation since it is proportional to R A , as compared to the 
R 2 A dependence from the first term due to LPM effect. 
We have to keep the first term for heavy quark propa- 
gation since the second term will have a similar nuclear 
dependence when the mass dependence of the gluon for- 
mation time is important. 

Since xl/xa ~ xbM 2 /xaQ 2 , there are two distinct 
limiting behaviors of the energy loss for different values of 
x B /Q 2 relative to x A /M 2 . When x B /Q 2 > x A /M 2 for 
small quark energy (large xb) or small Q 2 , the formation 
time of gluon radiation off a heavy quark is always smaller 
than the nuclear size. In this case, 1 — exp (— x 2 L /x 2 A ) ~ 1, 
so that there is no destructive LPM interference. The in- 
tegral in Eq. (14) is independent of Ra, and the heavy 
quark energy loss 
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is linear in nuclear size Ra- In the opposite limit, 
xb/Q 2 xa/M 2 , for large quark energy (small xb) or 
large Q 2 , the quark mass becomes negligible. The gluon 
formation time could still be much larger than the nu- 
clear size. The LPM suppression factor 1 — exp (— xj^/x^) 
will limit the available phase space for gluon radia- 
tion. The integral in Eq. (14) will be proportional to 
/ dxi[l — exp(— x\ /x 2 A )]/x 2 L ~ 1/xa- The heavy quark 
energy loss 



C A 



Ca 2 x B 
N c x\Q 2 
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now has a quadratic dependence on the nuclear size simi- 
lar to the light quark energy loss. Shown in Fig. 1 are the 



numerical results of the R A dependence of charm quark 
energy loss, rescaled by C(Q 2 )C'Aa 2 (Q 2 )/Nc, for differ- 
ent values ofxB and Q 2 . One can clearly see that the Ra 
dependence is quadratic for large values of Q 2 or small 
xb- The dependence becomes almost linear for small Q 2 
or large xb- The charm quark mass is set at M = 1.5 
GcV in the numerical calculation. 
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FIG. 1. The nuclear size, Ra, dependence of charm quark 
energy loss for different values of Q and xb- 
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FIG. 2. The Q 2 dependence of the ratio between charm 
quark and light quark energy loss in a large nucleus. 

Another mass effect on the induced gluon radiation is 
the "dead-cone" phenomenon [13] which suppresses the 
the small angle gluon radiation. Such a "dead-cone" ef- 
fect is manifested in Eq. (3) for the induced gluon spectra 
from a heavy quark which is suppressed by a factor 



fQ/q 
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_£ 2 T + z 2 M 2 _ 





(17) 



relative to that of a light quark for small angle radiation. 
Here 9 = M/q~ and 6 = £ T /q~z. This will lead to 
a reduced radiative energy loss of a heavy quark, amid 
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other mass dependence as contained in C{(z, M 2 ) in 
Eqs. (6)-(8). Setting M = in Eq. (14), we recover 
the energy loss for light quarks as in our previous study 
[18]. To illustrate the mass suppression of radiative en- 
ergy loss imposed by the "dead-cone" , we plot the ratio 
(Azf)(xB,Q 2 )/(Az<i)(x B ,Q 2 ) of charm quark and light 
quark energy loss as functions of Q 2 and x B in Figs. 2 
and 3. 
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FIG. 3. The xb dependence of the ratio between charm 
quark and light quark energy loss in a large nucleus. 



Apparently the heavy quark energy loss induced by 
gluon radiation is significantly suppressed as compared 
to a light quark when the momentum scale Q or the 
quark initial energy q~ is not too large as compared to 
the quark mass. Only in the limit M <C Q, q~ , is the 
mass effect negligible. Then the energy loss approaches 
that of a light quark. 

In summary, we have calculated medium modification 
of fragmentation and energy loss of heavy quarks in DIS 
in the twist expansion approach. We demonstrated that 
heavy quark mass not only suppresses small angle gluon 
radiation due to the "dead-cone" effect but also reduces 
the gluon formation time. This leads to a reduced radia- 
tive energy loss as well as a different medium size depen- 
dence (close to linear) , as compared to a light quark when 
the quark energy and the momentum scale Q are of the 
same order of magnitude as the quark mass. The result 
approaches that for a light quark when the quark mass is 
negligible as compared to the quark energy and the mo- 
mentum scale Q. Similar to the case of light quark prop- 
agation, the result can be easily extended to a hot and 
dense medium, which will have practical consequences 
for heavy quark production and suppression in heavy ion 
collisions. 
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